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The 2D screening is investigated in a simple single band square tight-binding model which qualita- 
tively resembles the known electronic structure in high temperature superconductors. The Coulomb 
kernel for the two particle Bethe-Salpeter equation in the single loop (RPA) approximation for the 
polarization can be evaluated in a strong tight binding limit. The results indicate an intense screen- 
ing of the Coulomb repulsion between the particles, which becomes stronger and anisotropic when 
the Fermi level approach half filling (or equivalently, when the Fermi surface approach the Van Hove 
singularities) and rapidly decreases away it. The effect is also more pronounced for quasi-momenta 
regions near the corners of the Brillouin cell, which correspond to dual spatial distances of the or- 
der few unit cells. Therefore, a possible mechanism is identified which could explain the existence 
of extremely small Cooper pairs in these materials, as bounded anisotropic composites joined by 
residual super-exchange or phonon interactions. 
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I. INTRODUCTION 

The small dimension of the Cooper pairs in the High 
Temperature (HTc) superconductors remains being an 
unclear point in the physics of these materialsiS^. In our 
view the main point to understand is the way in which 
the strong Coulomb repulsion which should be present 
between the two particles at distances of the order 15 
A (the estimated size of the pairs), should be somehow 
compensated for allowing the particles to be bounded. 
The cancellation in the old superconducting materials 
is carried out by the normal metallic screening at the 
Cooper pair sizes in these materials which are of the order 
of hundreds of A. It allows the weak phonon attraction 
to realize the binding. 

Considering this situation one is led to the idea that 
there should exist special characteristics of the HTc 
materials allowing the mentioned compensation of the 
Coulomb repulsion. One of the most salient features of 
the materials showing high Tc superconductivity is the 
presence of lattice planes formed by Cu and O atoms 
in a nearly square arrangement. Therefore, it can be 
naturally suspected that the 2D character of the elec- 
tron dynamics of those planes could be closely connected 
with the screening of the Coulomb repulsion at short dis- 
tances. Following this idea, in the present work we inves- 
tigate the screening of the Coulomb potential in a simple 
tight binding model chosen to qualitatively resemble the 
valence electron properties of the high Tc materials. In 
order to allow the arrival to analytical results, the model 
is constructed as simple as possible and purely two di- 
mensional. 

The main aim will be to obtain an estimation of the 
dielectric function and the kernel of the Bethe-Salpeter 
bound state equation associated to the Coulomb inter- 
action. The model will be a simple tight binding one 



designed, roughly speaking, to match the dispersion re- 
lation of the half filled band crossed by the Fermi energy 
in the high Tc materials. The energetic width of this 
band is extracted from the work Some additional 

simplifying assumptions also help to reduce the mathe- 
matical complications in the estimation of the dielectric 
function and the kernel. However, we think that they do 
not sacrifice the qualitative physical essence of the dis- 
cussion. It should be stressed that the justification for 
starting from this simple model, comes from the need 
of reducing the technical complications created by the 
lack of full translation invariance of the physical problem. 
This fact cancels the simple dependence of the dielectric 
response quantities of homogeneous systems on the differ- 
ence of the spatial arguments. This circumstance led us 
to search for analytical results in the non-interacting part 
of the electronic Hamiltonian allowing to better tackling 
the complications in evaluating the dielectric response. 
The calculations are done following the method devel- 
oped by Hanke and Sham for the calculation of the di- 
electric response properties of crystalline system a 5 i 6 i 7 . 

In our view the obtained results for the dielectric func- 
tion and the interaction kernel have interesting phys- 
ical implications. They indicate that, whenever other 
higher contributions do not balance out the effect of the 
here presented ones, the amount of screening produced 
by the planar electrons can be intense when the system is 
near half filling. Thus, the possibility is suggested of the 
occurrence of a strong screening of the Coulomb repul- 
sion in the real materials for small doping. Being close 
half filling, the evaluated dielectric functions can reduce 
the value of few eV of the bare Coulomb potential at 
distances of few periods of the spatial lattice, to values 
of the order of 0.1 eV. Therefore, a possibility is indi- 
cated for other weak binding mechanisms to furnish the 
necessary attractive forces for the pair creation. Cur- 
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rent ideas consider that a combination of diverse forces 
(super-exchange, electron-phonon, excitonic, polaronic, 
etc.) could play a role in the formation of the Cooper 
pairs. However, experiments have still not allowed to fix 
the nature of the acting mechanism, or at least, there is 
no consensus about the relevance of any particular onei^. 
It also follows that the screening effects rapidly decrease 
when the density of holes created on the half filled ground 
state increases. This property in competition with the 
natural need for a non vanishing hole density for the su- 
perconductor condensate to exist, could furnish an expla- 
nation for the first growing and after decreasing value of 
the critical temperature as a function of the hole density. 

Up to our knowledge, the idea that the Coulomb poten- 
tial in a 2D crystal system can contribute to produce, not 
only screening but even binding between pairs, was firstly 
advanced by D. MattiA Also various authors have been 
also considering the effects of screening in superconduc- 
tivity along various directions of thinking^iiiiiSiiiiiiii. 

A special point to be noticed is the fact that the here 
considered RPA approximation, in general grounds, is 
more suitable in situations in which the kinetic energy is 
greater than the interaction one— i. This means that the 
limit of validity of the here evaluated quantities should 
be closely examined, since the interactions are important 
in the high Tc materials. However, it can be stressed 
that the applicability of the RPA in Hubbard like mod- 
els in the metallic region up to near the metal-insulator 
transition has been argued 1 ^. This question is expected 
to be more closely addressed elsewhere. It should be 
also underlined that the results presented here support 
the also strong screening effects, and even over screen- 
ing, of the on-site Coulomb coefficients in single and two 
bands Hubbard models reported in [j| . The study of the 
possible connections between these results with the ones 
presented in this work seems to be worth considering in 
further studies. 

The particular motivation for realization of the present 
work was created by the already mentioned observations 
about the characteristics of the HTc superconductive sys- 
tems: a) The remarkable smallness of the Cooper pairs 
(of the order of 15 A), b) The foreseeable need of a mech- 
anism being effective in screening the Coulomb repulsion, 
which is of the order of few eV at distances of the size 
of the Cooper pairs. The existence of such a mechanism 
then could allow the short distance attraction of other 
participating forces, like the super-exchange interactions 
by example, to bind the pairsii. 

The structure of the work will be as follows: In Section 
2, the tight-binding model for the ab ceramic CuO planes 
will be defined. Simple Gaussian orbitals are introduced 
from the start, allowing (following Wannier 18 ) to find ex- 
plicit expressions for the Bloch functions. Then, the tight 
binding dispersion relation e(k) for the model is matched 
to the one shown by the existing half filled band in the 
HTSC materials 4 . In Section 3 the fermion free prop- 
agator of the model is employed to obtain formulae for 
the proper polarization and its Fourier transform in the 



ladder (RPA) approximation. Further, in Section 4, em- 
ploying the results for the proper polarization, a formula 
is also derived for the screened Coulomb potential. All 
these expressions show the complicate dependence on two 
arguments which follows from the lack of translation in- 
variance of the system. Then, Section 5 starts the appli- 
cation of the tight binding approximation to simplify the 
formula for the screened Coulomb potential. The final 
Section 6 is devoted to the calculation of close expressions 
for the dielectric function and the Coulomb screened ker- 
nel of the bound state Bethe-Salpeter equation. These 
quantities arise as functions of only one argument: the 
conserved transmitted reduced quasi momentum in the 
first Brillouin zone. This is the main technical result of 
work, since the breaking of the translational invariance 
of the problem obstacle the finding of such formulae in 
general. The simplifications were a direct consequence 
of the tight binding approximation and allowed to plot 
these quantities in the first Brillouin zone for reasonable 
values of the parameter defining the overlapping in the 
tight-binding model. The results indicate the possibility 
for the existence of a strong screening of the Coulomb 
interaction for conditions resembling the ones present in 
the real materials and they are commented in this ending 
section. 



II. THE SINGLE BAND TIGHT-BINDING 
MODEL 

This section introduces the simple model for the CuO 
planes in the superconducting ceramics. As the super- 
conductive properties along the CuO layers are signif- 
icantly stronger that in the perpendicular c axis, the 
movement will be considered as two dimensional. Then, 
it will be supposed that the electrons move in a square 2D 
lattice of unit cell size p — 3.82 A. This value of the unit 
cell parameter of the CuO lattice was borough from the 
workiS. In order to simplify the screening calculations 
in the considered non translational invariant problem, 
we will also adopt a picture motivated in the one band 
ij-model, in which simple Gaussian orbitals are centered 
only on the Cu sites. It is expected that the results for 
the screening effects on the Coulomb potential that are 
obtained in the work could be good indicators of what 
occurs in the real systems. 

The Gaussian orbitals staying on the points of the 
square lattice (See the figure ^) will have the form 



y>(x) = 



1 



N at = Tta 2 



(1) 



where r = |x| is the distance to the origin which will 
be taken as a fixed point of the lattice and a is a pos- 
itive constant. From now on, all the vectors will be 
represented in boldface characters. The function <p(x.) is 
normalized in the whole plane of the crystal. The value 
of a measures the width of the spatial region where the 
orbital have appreciable values. 
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FIG. 1: Picture illustrating the square lattice in which the 
Gaussian orbitals are centered. The unit cell size is p — 3.82A 
and the wave functions are denned on the plane. This sim- 
plification allows to obtain analytic results for some of the 
screening response quantities. The mean potential created by 
the ions and atomic cores is assumed to furnish the experi- 
mentally observed band width of the half filled band crossing 
the fermi level. This is another simplifying assumption done 
in the work. 

Let us now consider the determination of the ana- 
lytic form of the orthonormalized set of Bloch wave 
functionsiS* 2 ^, generated by the orbitals 0) when dis- 
placed on to all the points of the considered squared 
lattice. The valence electrons of this band will be as- 
sumed to be loosely bounded to the atoms laying on the 
lattice points. Then, they can propagate along the crys- 
tal feeling a periodic potential created by the ions and the 
rest of the electrons being more bounded to the atoms. 
These more localized electrons will be assumed to pertain 
to the filled energy bands of the whole crystal. Then, 
we will search for the Bloch functions associated to the 
valence electrons in the form 

* k (x) = ^=^e ik Xx-R), (2) 

R 

where R is the position vector of a general point of the 
square lattice. As usual, it follows that the functions 
satisfy the Bloch conditions 

* k (x + R')-e lkR '* k (x). 



The proposed Bloch orbitals will approximately de- 
scribe the valence electrons in the crystal receiving the 
action of the periodic potential of the atoms and core 
electrons. As stated before, the rest of the electrons are 
more deeply bounded to the atoms and their effect will 
be assumed to produce produce a kind of mean field 
influence on the valence electrons. 

The just defined functions become orthogonal among 
them due to the very same construction, since they are 
eigenfunctions of the translations in the lattice with dif- 
ferent eigenvalues. It rests only to define the normaliza- 
tion factor iV k in (J2J by imposing the additional validity 
of the orthonormality condition 

y"d 2 x* k ,(x)* k (x) = tf kik ,, (3) 

where (5 k , k < is the Kronecker delta in the indices k of 
the states. For this purpose let us substitute the sums 
defining the Bloch functions in and use the property 

5>* k - k ')- R = iV5 k , k ,, (4) 

R 

where N is the total number of points of the lattice and 
k = 2w(j^ni, j^Ti-i) for which L l5 L 2 are the spatial 
periods of the crystal in the two principal orthogonal 
directions, defined by the Born- Von Karman boundary 
conditions. After that, the normalization constant fol- 
lows in the form 



A k = N e~ ik R / <Pxp(x)<p(x - R), 

R ^ 

= N ]T e- ik R f d=V(x)T_ R ^(x) , (5) 

R J 

where Tr is the translation operator shifting the argu- 
ments of the functions in the vector R as Tr<^>(x) = 
¥>(x + R). 

The Bloch orbitals can be explicitly evaluated in terms 
of the Elliptic Theta functions. It can be noticed from the 
fact that their definition is given by a sum of Gaussian 
functions. Then, after substituting Q in J2J), it follows 



* k(x ) = * y e *-R _L e-^ 



1 



£< 



■ <"2°2-*2) +ifc2TO2a2 



e ^UlaT + ,qi)9 3 (— Tj + ^r- ,92), 



2a 2 i 



2a 2 i 



(6) 
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where it was defined R = aini + a 2 n 2 and x = a!x n i+ 
X2JI2 in which ni, 112 are unit vectors along the directions 
of the crystal axes. The precise definition of the Elliptic 
function employed was 

+00 

9 3 (u 7 q)= ?" 2 e 2m " • (7) 

n— — 00 

In an analogous way, the normalization factor can be 
calculated in the form 

N k = N9 3 (-^, qNl )6 3 (-^, qN2 ) . (8) 

The obtained Bloch orbitals satisfy the property 

*_ k ( X ) = n (x) , (9) 

which directly follows from substituting k by — k in © 
after taken into account the Elliptic function property: 
6* 3 (u,q)=8 3 (-u,q) . 

Let us impose now the tight binding approximation 
corresponding to the physical case in which the overlap- 
ping of the valence electron orbitals is small. Therefore, 
it will be considered that the wave functions are closely 
localized around the points of the lattice. Thus, the 
following approximation for the integral in JSJ) can be 
written 

J d 2 x^(x)f R (^(x) - 5 n ,o J d 3 V (x) = 6 n ,o ■ (10) 
But, substituting fTTijl in gives 

7v k = 7v]T e - 4k - R ^ R .o = 7v. (11) 

R 

That is, according to l|l(J|) . it follows 7V k m N and the 
Bloch functions take the simpler form: 

* k (x) = -^^e ikR ^(x-R) . 

R, 

Let us now consider the spectrum of the model, as- 
suming that the valence electron Hamiltonian H e will be 
the most relevant one for the description of the HTc su- 
perconductor properties. In other words, the action of 
the ions and the core electrons will be only considered 
as sources of mean external fields acting on the valence 
electrons. Therefore according to the Hartree Fock ap- 
proximation each valence electron will be moving in the 
periodic potential created by the interaction with the 
ions and the core atomic electrons. As the result of these 
assumptions, we will remain with a system having as its 
free Hamiltonian, the sum of the individual one particle 
Hamiltonians of the valence electrons 

H e = Tm + Uoi = ^ , (12) 
Hoi = Toi + Uoi , 



where is the Hamiltonian of a valence electron i. In 
addition, after also considering the tight binding approx- 
imation, it will be supposed that in the neighborhood of 
each point of the lattice, the free Hamiltonian of an elec- 
tron can be well approximated by an atomic Hamiltonian 
H at centered in the considered point. 

Therefore, under the above specifications the tight- 
binding procedure leads to the following dispersion re- 
lation for the valence electrons in the model under 
consideration^ 

e(k) = e F - 27 [cos{pk x ) + cos(pfc y )] , (13) 

where is the Fermi energy of the system and 7 defines 
the energetic width of the band as 87. 

Thus, the electrons will lay in states described by 
the wave functions ^(x) corresponding to energies 
given by 1|13|) where k = (k Xl k y ). Here the wave vec- 
tors k will chosen in the first Brillouin zone by conven- 
tion, since k'= k + Q correspond to physically equivalent 
states. With the help of (|13f) . it can be checked that 
e(k + Q) = e(k) as the equivalence implies. 

For the definition of our qualitative model, we simply 
matched the square lattice parameter to the experimental 
value. Also, it was directly assumed that the one particle 
potentials created by the ions and the atomic kernels are 
the necessary ones for fixing the width of the band in 
(|13fl to the estimated value. The lattice parameter was 
the one corresponding to the ab planes in the ceramic 
material YBaCuO. That is, p = a = b = 3.82 A as re- 
ported ini^. The energy width of the half filled electronic 
band in the material, as defined by the energy difference 
87 between the top and the bottom of the band, was es- 
timated here from the data in Ref— to be 3.7 eV. The 
relation (|13f) also well reproduce the form of the band 
given in*. 

Another assumption chosen for the sake of simplifi- 
cation is that the dispersion relation of the half filled 
band is solely determined by the ionic plus atomic ker- 
nel potential. It considers that the interaction between 
the valence electron does not strongly contribute to the 
band width, which is not necessarily valid Therefore, 
we expect in the further extension of the present work 
to perform a self consistent derivation of the band width 
starting from a more basic interacting Hamiltonian. 

Now, let us precise the form of the free propagator for 
the valence electrons to be employed in the further eval- 
uations. According to the above discussion, the electron 
states will have the form 

* kA (x) = * k (x)$ A , 

where Vt' k (x) are the Bloch wave functions and <&j = [1] 
y $| = [J] are the spinor states. In addition, as defined 

before, k = 2tt (j^, j^j where ni and n 2 are integers and 

Li,L 2 are the periods of the crystal along the principal 
directions. Then, the field operators can be written as 

$(x) = *ka(x) Qca , $ + (x) = ^ * kct (x)+ 3+ Q , 

ka ka 
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where a is the spin projection quantum number. There- 
fore, the free propagator can be evaluated from the usual 
formula^: 



iG a p{-xt,-x!t') 



(tto|r[$ ga (xf)$+3(x'tQ]|*I/o) 
(*o|*o> 



(14) 



where |Wo) is the non-interacting ground state of the 
model and T means the chronological operator. 

The free ground state to be considered will be one 
in which all the electron wave functions are filled up to 
some value of the Fermi energy ep . Since the con- 
sidered problem lacks the rotational invariance, there is 
not a unique Fermi wave vector kp associated to all the 
highest filled electron states. As usual, the field opera- 
tors can be rewritten in terms of the operators creating 
and destroying electrons above the Fermi energy or holes 



J 



below it, in the following way : 



s (x) = *kA(x)S kA + *kA(x)6" 



kA 



e(fc)>e F 



e(fe)<ej7 



(15) 



in which the usual definitions are given of the creation of 
holes operators: 6l k ^ as being equal to the annihilation 
operator of electrons at momenta k for energies lower 
than Fermi one. The new annihilation operators satisfy 

Ok/3 1 $o) = and fe-k/jl'&o) = , since there are no 
electrons over, nor holes below the Fermi level. Then, 
the expression for the valence electron free propagator 
can be evaluated following the usual steps in the standard 



k 



dui * kQ (x)* k/3 (x')e 



-iuj(t—t f ) 



2-7T hw — e(k) + ii].sgn(€(k) — cf) 



(16) 



III. PROPER POLARIZATION 



This Section will consider the evaluation of a formula 
for the polarization function. Further, these expression 
will be employed to evaluate the screened Coulomb po- 
tential and the kernel of the Bethe Salpeter bound state 
equation. The lowest order contribution to the proper 
polarization is given by the one loop term formed by 
two electron propagators. The associated Feynman dia- 
gram is illustrated in the figure |21 as the amputated loop 
and its analytic expression can be written as follows 



J 



P{x,x')xx>w> = G%(x,x') G%(x',x) , (17) 
= G°(x,x') GV,z)<VW , 
= P(x,x') <Wa' , (18) 

where for writing it, the formerly defined compact nota- 
tion x = (x, t) has been employed. 

After substituting (|16fl it follows 



P(x,x') = h 2 



duj f°° dtu' x - x - 

2^ / 9^2^L' 



^')(«-*')^ k ( x )^*(^)^ k ,( x ^^( x ) 
[hu -?(k)] [fiw' -?(k')] 



(19) 



r 



where it has been defined 



e(k) = e(k) - i rj sgn(e(k) - e F ) 



(20) 



As it can be observed, the polarization has explicit de- 
pendence on two spatial arguments. This is a natural 
outcome since the lattice potential breaks the full trans- 
lational symmetry. Therefore, the model under consider- 
ation does not have a continuous translational invariance 
and the technical evaluation of the dielectric properties 
turn out to be more cumbersome than for the homoge- 



neous systems. On another hand, since the Hamiltonian 
of the system is time invariant the time dependence of the 
polarization will be only through the time differences. 



Thus, let us determine the Fourier transform of 
P{x,x') over the two spatial arguments and the time. 
The Fourier variables will be associated according to the 
rules x — > q , x' — * q' and (t — t') — * ui e , and 
the convention to be used for the Fourier transforms of 
P(x, x', t — t') is explicitly defined by the relation: 
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FIG. 2: The first correction to the screened Coulomb po- 
tential. The relative difficulty in the present evaluation of 
this loop is produced by the lack of translational invariance 
of the lattice. In the momentum representation the impulse 
entering in the input line is in general different than the one 
in the outgoing line. However, the reduced translational in- 
variance of the problem will assure the conservation of the 
quasi-momentum as reduced to the first Brilouin zone. 



^(q,q',< 



= d 2 



d(t ~ t')e 



= (*- t ')F(x,x',<-i') 



Substituting P(x, x',i — t') by its expression (|19|l . con- Dirac delta definition 8(p) = ^- dx e l p ' x , results 
veniently reordering the operations and employing the in: 



J 



P(q,q',c e ) =^ 2 EE{/ ^ x J rf 2 x' e ^ x ^' x 'vl/ k (x)^(x')* k '(x')^(x)| 



k k 

00 dw 



_ x 2tt [hto-7(k)][h(Lo-LO-)-I(k<)Y 



(21) 



Let us decompose an arbitrary vector of the plane as 
the superposition of a translation vector in the di- 
rect lattice and a vector z contained in the elementary 
Wigner-Seitz cell. Then, for any vectors x and x' is 
possible to write 

x = z+[x] = z + R, 
x = z +|xj=z + R . 

The squared brackets [x] could be seen as an "integer 
part" of the considered coordinates in the plane. That 
is [x] = R where R is a lattice vector and z G W in 
which W is the elementary cell of the crystal. After that, 
the above expression for P(q, q', ui e ) can be simplified by 
substituting the integrations through all the plane over 
the variables x or x' by integrals over the elementary cells 



as: 

[d 2 xF(x)=[ d 2 zVr(z + R). (22) 
J Jw R 

After defining the quantity 

7= J d 2 xj d 2 x'e- lq x - 4q ' x '* k (x)**(x')*k'(x')*k'(x) , 

with the help of the Bloch functions property 

*k(x) - * k (z + R) = e ik R * k (z) , (23) 
and applying (|22|) . it can be written 
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I=( d 2 z ( dV* k (z)*J (z')tf k ' (z')*J, (z) cxp [-iq.z - iq'.z'] 
2 exp [i(k - k' - q).R] ^ cxp [-»(k - k' + .q')-R-'] ■ 



R 



R' 



r 



The following relation allow to simplify the last 
equation 



R 



,i(k-k'-q).R = N S (Q) ( k _ k ' _ q ) ; (24) 



where 6™ is the Kronecker delta in the space of the wave 



vectors k, defined in the first Brillouin zone. The super- 
index means that it is periodically extended outside the 
zone. That means 5^ (k - k' - q, 0) = 5^ (k - 
k' + q' + Q, 0), with Q being an arbitrary vector of the 
reciprocal lattice. Then, making use of (|24|) . the integral 
I reduces to 



I= jv 2 j( t W(k-k'-q,0)^ t »(k-k' + q',0) f d 2 z, f d 2 z'* k (z)^(z') 

Jw Jw 

*k'(z')*lk'( z ) ex P [-*q- z - iq'-z'] . 



r 



Substituting this relation in the polarization expression 
(j?Hl , grouping the terms having dependence on k, k' and 



after performing the sum over k', the following formula 
arises for the polarization loop 



P( q ,q^)=ft 2 A/ 2 -^^(q + q',0) r ^ I d 2 z / dVe-^-l'*' 

{^r J-oo 2tt j w j w 



(27T) 

J d 2 k* k (z)vl/ k (z')*k- q (z')* k - q (z) 



[Huj - e(k)] [H(u - u e ) - e(k - q)] 



(25) 



After performing the frequency integral over w, it also follows 



P(q,q',w e ) = (_^)JV 2 5(q + q') f d 2 z [ d 2 z' exp [-iq.z - iq'z'] / d 2 k* k (z)* k (z' 

?(e J? -e(k))^(e(k-q)-e F ) 



* k _ q (z')* k _ q (z) 



g(£ F - e (k-q))g(e(k)-ef) 

[Sw e + e(k - q) - e(k) + irj\ [hu e + e(k - q) - e(k) - irj\ 



(26) 



I 

where the periodic delta function <5(q) = jT-^^(q) na s e(k) = e(k + Q), implies e k , q +Q = £k,q), using the wave 
been defined. The limited lattice translational invariance property 5'k+Q(z) = ^k(z) and substituting q = q+ Q 
of the problem can be made explicit. Defining the the and q' = q' + Q' in 126( 1. the following alternative ex- 
quantity e(k — q) — e(k) = ek, q (for which the prop- pression can be written 
erties of the energy spectrum 1(13(1 e(— k) = e(k) and 
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P(q+Q,-q'-Q',w e ) = (--)jV 2 £(q-q') f d 2 z[ dVexp [-iq.z + iq'.z'] 

% Jw Jw 

exp HQ.z + zQ'.z'] J d 2 k* k (z)^(z')* k -q(z')*k-q(z) 

(e F - e(k - q))0(e(k) - e F ) 0(e F - e(k))0(e(k - q) - e F ) 



[huj e + e kyil + irj\ [hw e + e feiq - 

= ( 5(q-q')P(q,Q,Q> e ) 



IT]' 



(27) 
(28) 



The lattice invariance of the problem, is reflected in 
this relation by the delta function evaluated in the differ- 
ence between the reduced momenta. Thus, the limited 
translation invariance is able to assure the conservation of 
the reduced momentum in the loops, but it can't avoid 
the changes over the entering and outgoing reciprocal 
lattice vectors Q. 



The property H28fl will be of help for the simplification 
of the various terms in what follows. For its use is help- 
ful to eliminate the products of Heaviside functions. Em- 
ploying the property 9{x) = 1 — 9(— x) in the expression 
for P(q, q',o; e ) defined in Q26)). the product of Heaviside 
functions can be transformed in sums in arriving to the 
alternative formula 



P(q,q',tj e )= (-^J N 2 5(c[ + q') j d 2 z J d 2 z' exp [-iq.z - iqV] J d 2 k# k (z)* k (z') 

(e F -e(k-q)) 6{e F - e(k)) 



* k _ q (z')* k _ q (z) 



[hw e + e feiq + irf\ [huj e + e k ;, 



ii}' 



(29) 



The conservation of the reduced quasi-momentum can 
again be made explicit by writing this relation in the 
form: 

P(q + Q, q' - Q', uj e ) - <5(q - q')P(q, Q, Q', UJ e ) , 



where q is the reduced wave vector, and it has been de- 
fined: 



P(q,Q,Q',a, e )= - 



i I .iw 



j d 2 z / d 2 z'exp[-iq.(z-z')]exp[-iQ.z + iQ'.z'] / d 2 k 
Jw Jw J 

[e F -e(k-q)) 6(e F - e(k)) 



* k (z)**(z')* k - q (z')* k _ q (z) 



[huj e + e fci q + irj\ [huj e + e k ^ - irf] 



(30) 



IV. SCREENED COULOMB INTERACTION. 



represented in the form 



Let now consider the evaluation of the Fourier trans- 
form of the screened Coulomb potential, associated to 
the partially filled band model for the valence electrons 
in the HTc superconductors. The bare potential will be 



U Q (x 1 ,x 2 ) = V (xi -x 2 )a(ii -t 2 ) , (31) 

c 2 



K)(xi - X 2 ) 



47re xi - x 2 



(32) 



where the compact notation x = (x,i) continues being 
used in what follows 22 . Then, the formula: 



47T£ |X — X'l 



d 2 q 



(2tt) 2 2s |q 



£ giq.(x-x') 
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defines the Fourier transform of Vb(xi — X2) as Let's evaluate now the ladder approximation series for 

the interaction potential U(x,y) a i3s.y between two va- 
y r s _ e2 lence electrons. The first contribution to this quantity 

2s o |q|' is represented as a Feynman diagram in the figure^ The 

result have the expression 

Because the full translation invariance, the Fourier trans- 
form over the two spatial indices and the temporal dif- 
ference of the main magnitudes being considered (polar- 
ization, screened potential, etc.) will be taken. These 
transformations will be written for any such of these 
quantities, let us say Q, in the notation 

(33) 



U{x,y) af}Sl = U (x,y) a p Sl + (r)(-l) d 3 xi / d x' 1 U (x, xi) a pxx>P{x\, nn<Uo{x' x , y)^>s 1 + 



(34) 



(35) 



where Uo(x,y) a psy is the bare Coulomb kernel. In this 
relation Uq and P indicate the functional kernels clearly 
defined by the Coulomb potential and the polarization 
function in the first term of the expansion. The bare 
potential has the explicit form 

U (x, y) a ps 1 = Vb(x - y)5 a p8s 1 8(t x - t y ). 



Taking the Fourier transform of P according to l|33|l 
and evaluating the integrals for a general contribution, 
the following expression for the Fourier transform of the 
polarization function arises 



tf(q,q\a>«) = (2i) 2 «'(q + q')V„(q) + (i)(-l)V„(q)P(q,q>')V„(-q') + ■ 

P(-q 1 ,q2^ e )Vo(-q2)...P(-q N ,q',LU e )V (-q') , (36) 



where S p defines the Dirac delta function for the whole 
plane S p (q) = 5(q) Sq-q' and 5(q) is the already defined 
periodic delta function. 



V. TIGHT-BINDING APPROXIMATION FOR 
BETHE-SALPETER TWO PARTICLE KERNEL 

Let us pass now to implement the tight binding ap- 
proximation for the evaluation of the screened Coulomb 



potential. Consider for this purpose the Bloch theorem 
property ^k(x) = cxp(ikx) itk(x), where the itk(x) are 
periodic functions in the lattice. Applying this relation 
to the wave functions an d substituting the explicit 
form of the Gaussian orbitals £[J , it follows 
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* k (x) = -i=^e' k 'Vx-R), 



R 



N 1 ,., „ 1 ( (X- R) 2 \ 

expukx) ■ > exp (zk.iR — x) —= exp — 

WNk^ WWt \ 2a? J 



(37) 



r 



Now, let us assume that the Gaussian orbitals have almost will not depend on the wave vectors k. This 

very small overlapping integrals for nearest neighbor can be observed by substituting the first exponential by 

sites, allowing to consider a « p. Then, it can be the unit in u^(x). Assuming this approximation, i|ll|) 
written 



itk(x) 



1 v-^ .., .„ „ / (x-R) 5 
2^ exp (ik.(R - x)) exp -- 



R 



2a 2 



But, the values of the functions «k(x) will be signifi- 
cant only for |x — R| a and thus the exponential 
exp (ik.(R — x)) as 1. Therefore, the functions itk(x) 



writes according to 



u k (x) w u(x) 




Substituting now in the polarization expression H30|) 
the relation u k (x) rs w(x), gives 



P(q,Q,Q',uj e ) = (--) N 2 [ d 2 z u(z) 2 cxp[-iQ.z] f d 2 z'u(z') 2 exp [iQ'..z'] 
V 1 J Jw Jw 

(e F -e(k-q)) 0(e F - e(k)) 



(39) 



[hiu e + e fe .q + irj\ [hhj e + e fei q - ir)' 
I 



It can be noticed how all the exponentials in k and q 
were cancelled. We consider this outcome as the main 
technical advance in this work. It allowed to obtain a 
close expression for the interaction kernel and greatly 
simplified the evaluations. Now, let us consider the 
quantity 



p(Q) =N d 2 z u(z) 2 exp [- jQ.z] 



(40) 



which basically is the coefficient of the Fourier series ex- 
pansion of the periodic function (u(z)) 2 . For its calcula- 
tion let us substitute l|38|l in l|40|l to obtain 



p(Q) = N f d 
Jw 



NNn 



EE 

R R' 



exp 



(z-R) 

2a 2 
I 



exp 



2a 2 



exp [— zQ.z] 



Now, the double sum over the R is converted in a 
single one since the overlapping is considered as being 
a very small quantity. Thus, the overlapping integrals 
of orbitals being centered in different sites can be disre- 
garded. Realizing an estimate of how small should be a 
with respect to p for the overlapping to be considered 
small, it follows that when - = 0.33 then the overlap- 
ping integral is of the order 10 % and when ^ = 0.289 



the integral is 5 % . In the case of a = 0.5 A for which 
- = 0.13 the integral take the value 4.6 x 10~ 7 . In this 

way, up to values near below a = 1 A the approximation 
employed seems to be a reasonable one. The p(Q) can 
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be explicitly evaluated as 

p(Q) = 



(z-Rr 



exp [— iQ.z] 



— / cPz^exp - 
^at Jw " y 

= p(Q) = tt~ / d 2 x cxp f-^j) exp [iQ.x] 

Nat Jw \ a J 



exp 



-(Qi 



in arriving to which the change of variables z' = — z was 
made by renaming newly z' as z, and the property l|22|) 
and exp(iQ.R) = 1 have been used. 



The function p(Q) is real and after substituting it in 
(1391) the polarization function can be written as follows 



P(q,Q,Q',,/) 



\ [fiw e + e fej q + iry] [fiu; e + e fciq - ir)'] 



(2nf 



-- ) p*(Q)n(q,c>(Q') 



(41) 



r 



where the quantity 



n(q) = 



d 2 k f g( Cf . -e(k-q)) fl(e F -e(k)) 
(2tt) 2 I + e fc , q + iv] [fuJ e + £fe,q - W. 



has been defined. 

Now, let's determine the simplifications induced in the 
screened potential by the tight binding approximation 
being imposed. For this purpose, let us perform the 
change of variables q'i = — qi and also employ l|22|) . but 
now applied to the momentum variables as 

f d 2 qF(q) = f d 2 q^F(q + Q), 



where q = q+ Q,q is a reduced wave vectors and Q is 
an arbitrary reciprocal lattice vector. After considering 
that any of the polarization functions appearing has a 
Dirac delta in the difference between the entering and 
outgoing reduced momenta it follows 



C/(q,q> e ) = (2^) 2 <5(q-f-q')^o(q) + (2^) 2 ( 5(q+q')Vo(q+Q)p*(Q)n(q,^)p(Q')^o(q-Q^ 



53 .(27r) 2 ^q+q')y (q+Q)p*(Q)n(q,c i ; e ) 



N=l 



53/ 9 *(Q 1 )Vb(q+Q 1 )p(Qi) 

Qi 



n(q> e )x 



5>*(Qa)^»(q+Q 2 )p(Q2) 

Q 2 



.n(q,^ e ) 



X]/°*(QN)Vb(q+QN)p(QN) 

Qn 



n(q, W e )p(Q')K,(q-Q') 



(2^) 2 <5(q + q')Fo(q) + (2Tr) 2 6(q+q')V (q+Q)p%Qmq,u e )p(Q>)V (q-Q')- 

°° r in 
53 (2^) 2 <5(q+q')K ) (q+Q)p*(Q)E |n(q,W e )y(q)J n(q,u,>(Q')F (q-Q') , 



iV=l 



N=l 



where q and q' are the reduced quasi-momenta asso- 
ciated to q and q' and the sums over the Qi were 
accommodated in order to absorb all of them in the 
quantity that in what follows will be identified as the 



bare Coulomb kernel of the Bethe Salpeter bound state 
equation 

V(q) = $>*(Q)V (q+Q)p(Q). (42) 
Q 
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It is possible to incorporate in the second term in the 
last line of the last expression for U, in the appearing 



J 



summation to write 



U(q + Q,q' + Q', ^) = (2tt) 2 % + q')Vb(q) + (27r) 2 5(q+q')^o(q+Q)p*(Q) 

n(q,u, e )p*(Q)F (q-Q'MQ') 



N=0 



(2 7 r) 2 ,5(q+q')K ) (q) { <5q,-q' 



1 - n(q, w e )y(q) 



(43) 



where, also, the geometrical series was formally summed 



n(q)V(q) 



< 



over. This series is only convergent when 

1. As it will be the case, this condition will not be satis- 
fied in the whole Brillouin zone (See section 6). However, 
in our problem there will exist domains of the variable q 



over which 



II(q> e )V(q) 



< 1 and for them the summa- 
tion relation is obeyed. Therefore, after assuming that 
there is no change in the analytical expression for the 
considered quantity when the momenta is varied, the va- 
lidity of the resulting formula is taken as the analytical 
extension of the values in the convergence region. In sup- 
port of the above interpretation is also the fact that in the 
linear response theory 2 ^, the general posing of the prob- 
lem validates the employed formal summation formula. 
The equation Q43[) also can be written in the form 

U(q + Q, q' + Q', co e ) = (27r) 2 <5(q+q')^o(q)x(q, Q, Q'-^) , 

(44) 

just defining the quantity \- 



VI. BETHE-SALPETER COULOMB KERNEL 
AND DIELECTRIC FUNCTION 



Up to now, we have been considering the polarization 
and effective potential kernels of a non translational in- 
variant system. Therefore, the dependence of two spatial 
or momenta arguments makes their study more involved. 
However, it can be expected that the kernel of the two va- 
lence electrons (or holes) bound state problem, could be 
simplified by the conservation of the reduced momenta 
in the effective interaction. Thus, let us evaluate the 
kernel of the Bethe-Salpeter equation associated to the 
Coulomb interaction. This will be the relevant piece 
in the discussion of the effects of the screening to which 
this work is devoted. This kernel is determined by the 
matrix element of the screened Coulomb potential i|43|) 
describing the dispersion of two valence electrons. Its 
analytic expression is given by 



r(k,k',p, P ') 



d 3 3 



(45) 



The wave functions are taken in the interaction rep- 
resentation and the already defined convention is used: 
x = (x,i) and y = (y,t). They have the expressions 
— v I'k( x ) e ~ lLJkt ■ For this calculation the simpli- 
fying static limit u> e = will be assumed. That is, further 
we will only explore the screening in the static approxi- 
mation. 

In evaluating the matrix element the following steps 
were followed: a) The effective potential was substituted 
by its Fourier transform according to l|33[) . b) The trans- 
formation IL'l't was applied to the spatial x, y as well 



as the momentum integrals q, q', making use of the 



Bloch condition ^k(z + R) 



ik.R 



tion * k (z) 



Jk.R 



5'k(z) arid the rela- 



iik(z). c) The tight- binding ap- 



proximation (|38|) : that is, itk(z) « u(z) not depending of 
k, was implemented, d) The relation J2r e jk R = NS^.o 
was used, e) The standard formula for macroscopic 

crystals (5k, o = A S(k— 0) where A is area of the lat- 



tice, was employed, and finally: f) The following expres- 
sion for U was substituted 
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U(q + Q, q' + Q', 0) = (27r) 2 «5(q+q')^o(q)x(q, Q, Q, 0), 

n(q ! o)p*(Q)y (q-Q')p(Q') 



x(q,Q,Q',0) =*q,_q' 



l-n(q,0)V(q) 



(46) 
(47) 



as considered in the static limit u e = . 

After the above enumerated transformations the ma- 
trix clement (45) can be obtained in the following form: 



J 



(2tt) 3 N 2 

r(k, k', P , P ') = sfa - u k , + u p - uj p ,)S(k - k' + p p' y > A2 

E / d 2 z, it 2 (z)exp(iQ.z)E / d 2 zV(z') exp(iQ'.z') 

Q •* Q' J 

exp(-i(k - k' + p - p').z')Vb(k - k' + Q) X (q, Q, , 



r 



where x(q, Q, Q) = x(q, Q, Q', 0). 

Transforming the Dirac delta in a Kronecker one and 
using the representation l|40|l . after also considering that 



p(Q) is a real and even function, the previous expression 
is reduced to 



r(k, k', p, p') = S(p - p' + k - k', 0)<5(w k - + u p - cj p /) 

E E P(Q)p(Q')Vo(k - k' + Q) x (q, Q, Q) , 
Q Q' 



(2tt) 



<5(p - p' + k - k', 0)<S(c^ k - Wk/ + oj p - Wp ,)i^2r (q, k, k') 



r 



where the quantity To(q, k, k') is defined. 

Let us analyze the first term appearing in the def- 
inition of r (q, k, k') containing the Kronecker delta 
Jq^q', and which corresponds to the Coulomb poten- 
tial in the vacuum. For this contribution, it follows 

r«(q, k, k') = E p(Q)%(k - k' + Q) E P(Q')*Q,-Q' , 
Q Q' 

= Ep(Q)^o(k-k' + Q)p(Q) , 
Q 

= V(k - k') = y(q) . 



Performing a similar transformation for the second 
term appearing in i|47|) , which represents the modification 
of the Coulomb potential created by the 2D valence 
electrons, it follows: 
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IPS k. tf) . £ £ KQMQW(k - * + Q) HfiM^M93 

q Q' l-n(q)V(q) 

■ 2 n(q) 



(n(q)F(i 



(q)^(q) 



i 



n(q) 



1 - K(q)V(q) 
I 



Now, after summing the two contributions, the fol- 
lowing simple result arises for the Coulomb part of the 
kernel 



r (q, k, k') = if 3 (q, k, k') + I*"' (q, k, k' 



(2), 



n(q)F(q) 



n(q) 1 - n(q)F(q) ' 



^(q) 



1 - n(q)V(q) 



Henceforth, the dielectric function related with the 
screening of the bare potential has the formula 



e(q) = 1 - II(q)V(q) 



(48) 



The evaluation of e(q) was done for the overlapping 
parameter having value 

a = 0.95 A, 

and for a Fermi energy tp = —0.005 eV, being very close, 
but below the mid of the band. Thus, there is a very low 
density of holes in the system. The result of the calcula- 
tion plotted over the Brillouin cell is shown in the figure 
[3J As it can be observed the dielectric function shows 
high values for special zones in the Brillouin cell that run 
from more than 30 near the zone corners (like the point 
q = (p ' p)) down to less than 10 in other regions. Note 
that the zones near the diagonals correspond to rela- 
tive high values of the dielectric response. The fourth 
order symmetry of the result is also evidenced, although 
some numerical errors due to the presence of singularities 
in the evaluated integrals weakly break it. The peaks 
along the diagonals are only an artifact resulting from 
the evaluation for a finite number of points in the Bril- 
louin cell. The obtained results momentum dependence 
of the dielectric function becomes a consequence of the 
2D nature of the problem. Moreover, if the half filling 
condition is approached closely the dielectric constant 
increases drastically (up to values of the order of 2500 
were numerically evaluated by us). Therefore, it is clear 
that the Van-Hove singularities play a central role in the 
obtained effect. An additional and important outcome is 
that upon varying the Fermi energy away from half filling 




u i 

IMA" 1 ) 



FIG. 3: The dependence of the dielectric function on the 
quasi-momentum q in the Brillouin zone, for a value a = 0.95 
for the overlapping parameter. 



(e.F = 0), the values of the dielectric function along the 
diagonals rapidly decreases and the dependence becomes 
more isotropic. That is, the creation of holes reduces the 
screening of the Coulomb potential. This property could 
bring a natural explanation for the firstly rising and after 
decreasing of the critical temperatures in the HTc ma- 
terials. The picture could be as follows: first, when the 
number of hole density grows in the low density region, 
the strong screening state is established and the super- 
conductivity becomes stronger (Tc growing). However, 
as the hole density is increased even more the screening 
effect becomes weaker, and cancelling the enforcing ef- 
fect of the growing density by increasing the Coulomb 
repulsion. The rapid weakening effect given by the cal- 
culations done here support this last property. 

It should remarked that in this work we had not taken 
into account the dielectric constant of the medium Sb 
(due to the polarization of the ions and atomic cores). 
It seems that at distances of few unit cells this dielectric 
response could not be fully developed, but certainly its 
effects could be not so weak. Therefore, it seems rea- 
sonable to introduce a multiplicative factor A ranging 
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between and 1, describing the effective proportion 
of the dielectric response of the ionic and core medium, 
which could be acting at small distances of the order of 
the Cooper pair size. For large distances, the parameter 
should be equal to one and for very short ones should 
tend to become smaller than one. In general, it could 
expected that this effect can contribute even more to the 
screening of the Coulomb interaction obtained here. 




sidered approximation, the screening effect is noticeable. 
Considering that the bare Coulomb interaction produces 
a repulsive potential of few eV at distances of the order 
of the lattice constant, it follows that the screening effect 
could reduce these values down to ones of the order 0.1 
eV . But, at such small repulsion forces other mecha- 
nisms, such as super-exchange or strong phonon interac- 
tions could then create the necessary binding forces for 
the pairs to form. 

Finally it should be remarked that in order to ap- 
proach the present discussion to the real situation in the 
superconductor materials, it seems necessary to derive 
a similar picture but in a self-consistent approach. A 
concrete way for this program could proceed in the fol- 
lowing steps: a) To formulate the general Coulomb in- 
teracting problem of the valence electrons but retaining 
the Coulomb interaction exactly, b) Attempt to derive 
the here discussed tight binding model as a kind of HF 
approximation of the exact theory, c) After that, to de- 
rive and solve the Bethe-Salpeter equation for bounded 
pairs, which should be expected to end including the 
anti-ferromagnetic interactions of the tj-model as a pos- 
sible bounding mechanism^. The consideration of this 
program is in progress. 



FIG. 4: Plot of the screened and bare Coulomb contributions 
to the kernel of the Bethe-Salpeter equation, as functions of 
the wave vector q. Note the presence of an intense screening of 
the Coulomb repulsion in the considered here approximation. 

Further, in the figure 0] the results of the evaluation of 
the screened Coulomb kernel of the bound state equation 
for two electrons (or two holes) are plotted, in common 
with the values of the bare kernel. Again the overlap- 
ping parameter and the Fermi energy are a — 0.95 A, 
€p = —0.005 eV . As it can be observed, in the con- 
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